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ABSTRACT 


КО Ое ОШ AMC IIIE control systems is the subject 
of much emphasis at the present time. Test data collected 
ПКО С ИЕСИ resi conducted by the United States 
Army Material Systems Analysis Activity is examined using 
several different analysis procedures. The objective of the 
analysis is to determine the optimal time to employ in the 
ЕСЕТ ае ии тег ес а тапк тіге control lead prediction 
Booman л ас топа Objective is to quantify individual 
components of target error and to investigate the character 
EEG еггаге ав а Тіпстіоп от varying tracking rate 


filter times. 
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|. INTRODUCTION 


Currently much emphasis is being placed on the develop- 
ment of advanced tank fire control systems. Tanks incorpo- 
rating these advanced systems are expected to constitute a 
САКЕ КОЗА аА ИОТ the Allied force opposing the formidable 
Soviet threat present in Eastern Europe. 

This thesis presents an analysis of data collected dur- 
ing a tank fire control test conducted by the United States 
Army Material Systems Analysis Activity (AMSAA) in September 
1977. The specific area of analysis involves the investiga- 
tion of the optimal time to employ in the tracking rate fil- 
БП а rank fire control lead predicitTon algorithm. 
The analysis also had the objective of identifying individual 
components of tank gunnery error. An investigation was made 
of the relationship between these errors and the value 
uc £ ipe Tracking rare filter utilized in the fire 
٣٠٠٢٣٢٢٢٢ ٴً‎ 1. 

Section l| of This thesis motivates the need to pursue 
the development of advanced tank fire control systems and 
provides a historical development of improvements which have 
been made in fire control technology. An overview of recent 
experimentation which has been conducted to investigate The 
dynamics associated with engaging moving targets is included. 

The test purpose and test scenario is described in suffi- 


cient detail in Section lll to provide the reader with an 





understanding of how and why the test was accomplished. A 
Knowledge of the test scenario allows the reader to evaluate 
The analyses, and rhe analyses results, in light of the re- 
seu dcueomseimEosed Dy The scenario. 

The test data obtained from AMSAA was the raw instrumen- 
tation data as it was recorded during the conduct of the 
test. Before any analysis could be performed on the data it 
had to be expressed in standard units of measure. A descrip- 
Medea reduction process is included in Section IV. 

Initially four analysis procedures were used to investi- 
gaie ihe iest data. Procedures | and 2 dealt with analyses 
whose objectives were to establish a relationship between 
feaneazimush tracking error and varying tracking filter inter- 
vals. These intervals were called "delta time". Procedures 
5 and 4 investigated the relationship between the mean azi- 
muth tracking rate error and delta Time. 

The final phase of the thesis analysis used two models 
٠۰٠٠۰-٣٣۹٣۹3 target azimuth error as a function of delta 
Time. Through the use of these models an optimal delta time 
ines tracking Fate filter was identified. Additionally, 
Valuable insights were gained relative to the characteristics 


End dual components ot total target azimuth error. 





Il. BACKGROUND 


I STORTCAL DEVELOPMENT 

The events of the 1967 Arab-Israeli War once again served 
to remind the nations of the world who maintain standing ar- 
mies that mobile armored warfare can be, and often is the 
element within the combined arms team which can most effec- 
tively bring about success on the battlefield. Not since the 
early campaigns of World War 11 in which the German panzer 
divisions introduced Blitzkrieg, or lightning war, had the 
effectiveness of tank warfare been as devastatingly empha- 
sized. 

٣٠٠٠٠٠٦٣٦٠٦٠٦٦٦٦ Dasic tactics Of armored warfare as advanced 
by Liddell Hart, Heinz Guderian, and Erwin Rommel have 
changed little since inception, the same is not true of The 
fighting machines employed in the armored conflict. The mod- 
ern battle tank bears little resemblance to its forerunners 
of World War Il. Dramatic technological advances have been 
made in areas such as metallurgy which has vastly enhanced 
the effectiveness of armor plating. Development of more re- 
liable and durable vehicle power plants and suspension sys- 
КЕШЕ е he medern tank mobility and durability character- 
istics which would have been difficult to envision three de- 
cades ago. Similar improvements have occurred in the areas 
Serine power and fire control. Clearly all of these advance- 


ments serve To leave the modern tank with little in common 





with its forerunners except perhaps for the manner in which 
one goes about evaluating tank effectiveness. It has been 
said that of all battlefield targets which the tank is ex- 
шекеге епсасе,) The MOST difficult consists of hostile 
Tanks. As a consequence, the ability to "kill" enemy tanks 
has become the criterion of the tactical effectiveness of 
Tanks. The adoption of this criterion implies that tanks 
which are effective against enemy tanks are at least as ef- 
fective against targets other fhan tanks. While this may 
not be universally true, in general the ability to destroy 
other tanks is a satisfactory baseline measure of tank effec- 
٢٣ )٣ ة٠‎ 

Unlike tanks of the past, tanks on the modern battlefield 
must be able to successfully engage and kill moving enemy 
tanks in order to ensure their own survival. Although the 
ability to successfully neutralize moving enemy targets has 
always been a highly desirable capability to possess, its 
utter necessity is of relatively recent vintage. Technologi- 
tou co5neNghs in the area of turret stabilization now per- 
mit the modern tank to fire its main gun at the same time 
СИЕ ЕеБи А its target. In order to survive, The 
Tank being engaged must be able to neutralize its moving 
adversary. 

The key to the delivery of effective neutralizing fire 
against a moving target is a fire control system which effec- 
Tively assists the tank crew in the employment of the 
КЕШСЕ sedumamenhi. It would be difficult To locate an ex- 


perienced tank gunner who would not concede that one of The 
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more difficult aspects of tank gunnery is the precise esti- 
патоп of range to the target. To assist the gunner in range 
defermination, the United States Army began installing opti- 
cal range finders on its Tanks soon after the close of World 
Teen i. ЕОС па the installation of these devices no 
dramatic new advances were made in the area of range determi- 
nation until the recent advent of the laser range finder. 

The laser range finder's ability to very accurately determine 
range to target with little or no degradation by ambient light 
or weather conditions makes it an ideal component of any 
state-of-the-art fire control system. Range information is 
Produced by the laser range finder in the form of electrical 
impulses. These impulses become meaningful to the tank gun- 
ner only after they have been properly interpreted by another 
element of the fire control system, the electronic ballistics 
computer. Although the electronic ballistics computer serves 
the same basic function as the mechanical computer which itf 
replaced, its role in the fire control system has been great- 


ly expanded because of the order of magnitude increase in 


computing capability which it possesses over its predecessor. 
With input data from ancillary fire control system components, 
fhe electronic ballistics computer can compute target range, 


vehicle cant, crosswind, air temperature, powder temperature, 
air pressure, and tube wear parameters. Those parameters are 
in turn used by the computer to determine corrections in main 
gun azimuth and elevation to compensate for the effects of 

ШЕЛІ ЧГ рага ах, дип jump, tube droop, and crosswind 


on the trajectory of the projectile. Included in the azimuth 





СОЛИСТИ 15 The lead required to hit a target which is 


moving [2]. 


ER RECENT EXPERIMENTATION 
Within the past several years there has been a high pri- 
ority among members of the armor community and material devel- 


Spiers to determine the relationship, if any, between combat 


oie Survival on the battlefield and vehicular mobility 
and agility. Mobility as used within the present context re- 
fers to the ability of a vehicle to reduce its exposure time 


To a hostile force by moving at a high speed. Agility is de- 
fined as the ability of a vehicle to perform evasive maneuvers. 
It should be clear that any increased survivability which a 
combat vehicle gains because of its inherent mobility or 
Peery геси! т< from The fact that these characteristics have 
lessened the effectiveness of the fire being delivered by The 
hostile combat vehicle. Alternatively expressed, the diffi- 
Em encountered by the hostile vehicle in hitting its in- 
tended target may be increased because the effectiveness of 
its fire control system has in part been neutralized. Having 
Thus hypothesized the foregoing implicit, albeit subtle rela- 
tionship between degradation of fire control system effective- 
ness and survivability, it seems desirable to attempt to gain 
БООКЕ Уе іптогпатісп regarding their correlation. 

Experiments conducted in the Federal Republic of Germany 
iaa ys is Of data collected af The United States Armor Cen- 
ter have indicated that high mobility/agility may increase 


comba уепісіс survival on the battlefield. Specifically тһе 
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German experiment indicated that there was a high survivabi- 
MAD Ori Wien vehicular lateral accelerations of approxi- 
mately 0./g were attained. Two field tests recently con- 
ducted in the United States gathered gunner tracking perform- 
ance and hit probability data against maneuvering targets. 
These tests were the S-Tank Agility/Survivability (STAGS) 
Test and the Anti-Tank Missile Test (ATMT). Because of test 
Condi Ons and other constraining factors, neither of these 
Tests generated vehicular lateral accelerations above 0.4g; 
Mie Feroe, The data from These tests could not be used to 
substantiate earlier findings concerning vehicle survivabili- 
E 6 to lateral accelerations in the 0.7g range [3]. 
The most recent effort to gather field test data regard- 
ing high mobility/agility vehicular characteristics and their 
corresponding correlation to survivability has been under- 
Taken as a partial objective of the Armored Combat Vehicle 
Technology (ACVT) program. In order to accomplish this objec- 
tive, the United States Army Combat Developments Experimenta- 
tion Command (CDEC) was given a directive to conduct a High 
Mobility/Agility field experiment beginning in November 1977. 
The experiment was named the High Mobility/Agility, Phase IIA, 
Extended (HIMAG IIA EXTENDED) and was conducted at The CDEC 
Mies SITE ai Fori Hunter Liggett, CA. The HIMAG IIA 
test utilized a wheeled target vehicle which performed eva- 
Sive manuevers (designated sinusoidal wave patterns which 
generated lateral accelerations up to 0.7g) on an airfield 
Шоо at the field test site. During The کے‎ The tar- 


get vehicle was tracked by TOW and M60AI gunners who were 
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positioned at specified ranges and offsets relative to the 
vehicle maneuver path. Analysis of the gunner tracking data 
relative to the apparent motion of the test vehicle was in- 
ie need no Provide intiermation regarding the relationship be- 
tween vehicular mobility/agility and vehicle survivability. 
As indicated earlier, the testing was conducted to provide 
data from which a relationship (Cif any) between degradation 
of fire control effectiveness and target vehicle mobility/ 
agility can be established. 

In September 1977, the United States Army Material Systems 
Analysis Activity concluded a test at Aberdeen Proving Ground 
(APG), MD, whose objectives were very similar to those set 
grec ce HABMAG IA TesT. The testing at APG differed 
from that planned in the HIMAG 11А test in that the APG test 
was conducted in a highly instrumented laboratory environment 
and therefore could not be classified as a true field Test. 
E 2-0 )sweengueted in This thesis utilizes a portion of 
sx ucolheceted In the APG Test. A description of that 


test is presented in a subsequent section of this thesis. 
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۰٦٦ 2... DESCRIPTION 


As discussed in Section |l, the dynamics of the modern 
battlefield place a high premium on the modern tank's 
ability to successfully engage hostile moving targets. Future 
fielding of vehicles which possess increased mobility/agility 
Mor acteristics will make the task of hitting these highly 
maneuverable targets more difficult and will therefore inten- 
sify the need for more effective fire control systems. The 


Testing completed by AMSAA at APG in September 1977 repre- 


sents an important link in the development chain necessary 
to field a fire control system which will meet future chal- 
lenges. 


PURPOSE 

The APG test was entitled "Supplemental Fire Control 
Test, M60A3" and its purpose was to examine the implementa- 
ШЕШ Ке гүге control lead prediction algorithm, and To 
assess the system's capability to cope with moving targets 
traveling at varying speeds and directions of motion. An 
ЗИКА by definition, a rule for solving a certain 
type of problem. The problem to be solved in this instance 
Ca O 810110, To The maximum extent possible, a tank 
gunner in hitting hostile targets which are moving as pre- 
viously described. ln the absence of a fire control system 
Ме согрогатеѕ а lead prediction algorithm, the Tank 


Gunner is required to mentally compute the lead necessary to 
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ӨТІС ПО/ІИа targeT and having computed this lead, he must 
manually move the turret controls causing the lead to be in- 
۰۰۱۱۰۱٠۰۱1٦۰٦٦٦ enning azimuth. Using this method, the train- 
ing, experience, and aptitude of the individual tank gunner 
are key elements in whether or not he is successful in hit- 
mS Phe moving target. 

ШОО АРЫС tests were conducted using an M60AIE5 fire con- 
trol system. This fire control system assists the tank gun- 
ner by automatically computing the lead required to hit a 
moving target and by automatically moving the turret, thereby 
ieconporaring The computed lead into the firing azimuth. The 
lead computation is made by the fire control ballistics com- 
Mi ccosding TO The lead prediction algorithm. A typical 
moving target firing sequence for a tank equipped with an 
M60AIE3 fire control system would be the following. The tank 
Gunner places the center of the sight reticle pattern on the 
center of mass of the moving target and tracks in this manner 
tor a specified amount of time. After the gunner has main- 
Tained his track of the target vehicle for the minimum amount 
of time he engages a lead lock fire control switch which ini- 
Tiates the ballistic computer computation of the lead. Having 
ЕШШ Еле lead, The fire control system automatically in- 
Ее Ше ел Dy causing ne tank turret To rotate by an 
ЗЕЕ ато Ш Ч“ПІСП соггевропсв To the calculated lead. AT 
the same time that the turret is being rotated, the sight 
reticle pattern is moved by the amount of lead in the direc- 
tion opposite to that in which the turret is rotated.  Move- 


ment of the sight reticle pattern allows the gunner to 
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continue to track the target vehicle center of mass even 
гпошал in actuality the target vehicle is being led. The 
ОКЕ ageso ered bY a fire control system which incorpo- 
Киз Tama uromati¡c lead option can be summarized as follows. 
٣٣۷٦ا‎ ٥۸۷ ٣٠٠٠٠ "16" ۱97196٤6۴ required to compute the neces- 
sary target lead mentally. In essence, the human error in- 
pcdcs lead Calculation is eliminated. Second, the fire 
control system automatically moves the turret in azimuth by 
an amount corresponding to the computed lead. Since the gun- 
۰۰۰٠٠۱۱۱٦ ٣۷٣0٦٦٦٣ The responsibility of manually moving the 
Icon Trolls: ne is more able to concentrate on tracking 
The maneuvering target vehicle. Finally, the movement of 

the sight reticle pattern allows the gunner to continue to 
track the center of mass of the target vehicle versus having 
0+0 ۹۹۹ ی۶‎ (٥( the sight rericie pattern at a constant distance 
(actually at a constant angle which corresponds to the calcu- 
БО ела) fin from of The target vehicle. 

Based on the foregoing discussion, it seems clear that a 
fire control system which automatically computes and inserts 
target lead in moving target situations is a very desirable 
option to pursue. Given that premise, the next logical ques- 
tion to be answered is the following. What are the critical 
elements Of fire control lead prediction algorithm which 
when optimized, yield the best solution to the problem? Al- 
though there may be alternative solutions to this question, 
the analysis performed in this thesis deals with an investi- 
5 of the optimal characteristics of the tracking rate 


INE er. 
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ОШО СОГОТ OT ihe function performed by the tracking 
Tenens helptul in illustrating the vital part it 
Dee еа prediction algorithm. The lead prediction 
ВС От relies on the fact that when the gunner is tracking 
bie mMmOVing target through the fire control sight, the tra- 

s r [efor тле Turret corresponds to the angular rate of 
Senge от Тпе Target as viewed from the firing platform (the 
Tank). When the correct distance to the target is known (as 
determined by the laser range finder) and is coupled with 
other factors such as projectile type, powder temperature 
a UbDe wear, the ballistics computer can utilize this data 
To compute the projectile time of flight (TOF) from the gun 
to the target. The TOF multiplied by the mean angular tra- 
verse rate of the tank turret equates to the angular lead 
Redio nit the moving Target. From the relationship, 
Lead = TOF x Mean Angular Rate of Turret 
Mie learn that The correct specification of the angular 
ШО ШОШ се ©0 8 correct lead computation. 

By definition, rate represents change per unit of Time. 
The question arising in the foregoing problem centers around 
choosing the "unit of time" which will be optimal in terms 
SANZ ina the error Of The lead prediction algorithm. 
uu Tor те as specified in This context will be re- 
ferred to in the remainder of this thesis as "delta time". 
The problem can be illustrated by using examples which repre- 
sent extreme cases. Assume that The angular turret rate was 
computed based on the mean traverse rate of the turret over 


a delta timeof 10 seconds. Suppose that the target vehicle 
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was not moving for the first 5 seconds of the delta time per- 
٠٠ e Timo Ned at a mean rate of 3 mils per second thereafter. 
The mean angular rate of movement in this example is 1.5 mils 
per second; however, this is not representative of the move- 
emt sot the target in the critical period just before firing 
takes place. A similar extreme case can be visualized at the 
Opposite end of the spectrum whereby the mean angular rate is 
computed based on a delta time period which is too small fo 
5e realistically representative of Тһе true angular rate. IT 
seems apparent that the optimal delta Time parameter lies 
somewhere between large delta time values (10.0 seconds will 
be considered large) and those approaching zero as a limit. 
The purpose of this thesis is to present an analysis of the 
APG test data with an objective of isolating the range where- 


in optimal delta time values occur. 


ol) PROCEDURE‏ ہے 

The APG fire control test under consideration was composed 
of two phases, a steady-state phase and a transient phase. 
Each phase involved simulated firings by a stationary M6OAI 
ШАШКЕ е гт a MOOATES turret, at a simulated target 
which was moving in a horizontal plane. Because of the non- 
Manoa epe n Or the fest, if was accomplished in a building 
zm eOorespeciaily equipped for The conduct of fire control 
testing. The terminology "steady-state" refers to the motion 
patera Whieh he simulated target described. The Test was 
designed so that the evaluation of the fire control lead pre- 
diction algorithm would be based on the data collected during 


The steady-state test phase; whereas, other types of analyses 


19 





ШЕШЕ 56 accomplished using the transient test data. Because 


mussubwecr cor this Thesis involves an analysis pertaining to 


ШИЕ ЕС ОПТ О ead prediction algorithm, a full test de- 
scription will be presented for the steady-state test phase 
0۷ہ‎ 


Throughout the test, the simulated target was a laser- 
generated spot of light which was projected onto a large rect- 
angular screen positioned 31.5 meters in front of the center 
mass of The tank test vehicle. The motion of the simu- 
lated Target was controlled by projecting the output of the 
laser onto a mirror which in turn reflected the energy onto 
the target screen. The mirror used in the projection process 
was mounted on a drive mechanism so that the mirror motion 
could be controlled by a programmable mini-computer. The de- 
sired target motion was accomplished by programming the mini- 
Сре ит an appropriate function describing the motion. 
In the steady-state test phase the mini-computer was pro- 
grammed with a sine function resulting in a sinusoid profile 
being projected on the projection screen. As viewed through 
MI е О ОТ, mie back and forth motion of 
the projected laser spot represented the apparent motion which 
would be presented by a target vehicle moving at a constant 
ЕШ а сітепізг path. Simulation of varying target 
ranges and speeds was accomplished by varying the range and 
rate of the mirror angular movement. The steady-state test 
phase was conducted by running a sequence of tracking runs 
(simulated firings) at simulated ranges of 500 meters and 


1500 meters. Simulated target speeds were 50 kilometers per 
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hour (kph) at the 500 meter range and were 25 kph, 35.4 kph 
and 50 kph for the 1500 meter range. These range and speed 
combinations correspond to target vehicle maneuvers which 
wemrd generate lateral accelerations of 0.1/5g, 9.359 and 
ШЕП respectively. An illustration of the target simulation 


emebatrus and setup is depicted in Figure |. 


SIMULATED 
<— TARGET 
ج5‎ 








SCREEN 


N N 


COMPUTER 
Figure | 
The M60A| tank was emplaced with cant blocks under the 
miciieinack wich restilted in a 90.2 mil (5.07 degrees) cant. 


E unbese er the intentional cant was To cause the Target 


to have an apparent motion in elevation, as well as in 
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azimuth, even though the simulated target was being projected 
ОО ٦٦۹ ٦٦١٦۹۸٠۶ across the projection screen. Instru- 
mentation on the test vehicle included a gated television 
(TV) camera which recorded the view through the main gun fire 
Somo! sight and also a gated TV camera mounted on the gun 
muibemwhich was canted to correspond to the boresight of the 
main gun. Each gated TV camera was connected to a videotape 
recorder. Utilization of the gated TV instrumentation in con- 
ETON WIT a Target simulated by a high energy light source 
made it possibie to continuously record the position of the 
Be control sight and the position of the gun tube relative 
Ne POSTION Of target Throughout the test tracking runs. 
Data recorded in this manner were readily convertible to 
angular deviations and therefore became a fundamental portion 
of the data base for each of the test tracking runs. 

During the testing, the test system was operated by either 
Se avoetank gunners from The United States Army Armor Center, 
ПОИ По Kentucky. The Training and proficiency of these 
individuals was considered by the sponsoring agency at the 
КЕПӘР Genrer tO be representative of Typical Tank gunners 
found throughout the Army. Each gunner was allowed to gain 
familiarity with the test equipment and the test scenario 
ЕРЕ ОЛУ С ПЕЦ еа in a tracking run (being recorded) 
fer record, 

Within each recorded tracking run there were approximate- 
Рр осе firings. the gunners were instructed to 
track the simulated target during the entire duration of a 


given Tracking run. The test system was set-up so That The 
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test director could initiate a simulated firing sequence by 
s j i G a ready light located in the gunner's sight. When 
Enesgunner observed the illumination of the ready light he 
Mulli ied The following firing sequence. Tracking of the 
simulated target was maintained by the gunner by his slewing 
ИЕ Тиггет апа дап Cand correspondingly the fire control 
Signt) at the appropriate rate in azimuth and elevation. 

When the gunner made the judgement that his tracking motions 
matched the motions of the target he activated the fire con- 
trol lead lock enable switch. Activation of the lead lock 
enable switch initiated the computation of required lead and 
(15 subsequent automatic insertion into the firing azimuth. 
After the lead insertion and simultaneous movement of the 
sight reticle pattern, the gunner reacquired target tracking. 
ШИЕ ТІГІП sequence culminated ¡in the gunner's activation of 
Mening mechanism. throughout the firing sequence, test 
instrumentation recorded critical data relative to the command 
to fire, the activation of the lead lock enable switch and 
the gunner'!s initiation of a trigger pull. The commands to 
NEC ccm vcemhhessumnnempMwere given at various position loca- 
Bons or the Target on its sinusoid path so that a full spec- 
trum of its apparent motion would be represented in The re- 


GCG tracking run data. 
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IV. DATA REDUCTION METHODOLOGY 


0 INTRODUCTION 

lil sc cC i п сопа пе а detailed discussion of the data 
geduction methodology. It is provided for those readers not 
familiar with the character of the Supplemental Fire Control 
ET dala. Those persons not in this category may skip this 


Een without loss of continuity. 


B. RAW DATA TRANSFORMATION 

A characteristic of many highly instrumented testing pro- 
Geass Thal the dara collected during The conduct of the 
testing is very seldom useful for analysis purposes in its 
raw form. Most often, before any analysis can be conducted 
The raw test data must be converted into units which are 
meaningful in terms of the subject under investigation. The 
Esos applicable To the fire control test data which 
is analyzed in this thesis. 

The APG test consisted of a series of separate tracking 
Шет Edemeinarvyidual test tracking run was characterized by 
a designated run number, target range, target speed and gun- 
ner identification. Approximately 32 runs were conducted 
Шипа Іле steady-state Testing phase. Each tracking run 
Was roughly 4 minutes in duration. The analysis performed in 
this thesis used data from seven tracking runs. The remainder 
of the tracking runs were not analyzed because the time avail- 


able did not permit both the development of analysis method- 


24 





sy and analysis of all the test data. In light of The 
Time constraint, primary emphasis was placed on the develop- 
ment of analysis methodology. All of the runs selected for 
Aa Sis (гапв 166 through |42) were conducted at a simu- 
lated range of 1500 meters and with a simulated target speed 
of 50 kph. This range and speed combination corresponded to 
maximum target vehicle lateral accelerations of 0.7 g. As 
noted previously, there is much interest in obtaining data 
relative to gunner performance against targets capable of 
generating high lateral accelerations. The selected target 
range and speed combination was chosen on the basis of the 
foregoing consideration. Although only one range and speed 
Poeemeenvestigated in this thesis the analysis methodology 
employed is equally applicable to tracking runs conducted 

at different simulated speeds and distances. 

Throughout each test tracking run, data from the test 
instrumentation was collected and recorded on magnetic tape 
at an average rate of 106.55 samples per second. The compo- 
Sition of an individual magnetic tape data frame, or logical 
mecorg. included the following elements of information: Time 
the logical record was created; tracking/firing sequence 
event code; mirror drive position location; and azimuth/ele- 
vation position location relative to the simulated target. 
Or her Types of data were also included in the logical record; 
however, that data was not used in the present analysis and 
therefore will not be discussed. 

Шие селега)і approach taken with regard fo the data trans- 


formation was to write a Fortran computer program which would 
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transfer selected portions of the raw data onto a blank mag- 
петіс Tape while simultaneously accomplishing desired trans- 
formations on the raw data during the transcription process. 
The magnetic tape created as a result of this procedure con- 
tained the test data expressed in a meaningful form for ana- 
lysis purposes. Subsequent paragraphs will include a descrip- 
Tion of rhe character of the raw data and of the transforma- 
fos which were performed during the tape transcription pro- 
cess. 

As previously stated, logical records were created at an 
average rate of 106.55 records per second. The time component 
ell logical record was recorded in units of hours, min- 
utes and seconds to the nearest hundredth. Because of the 
manner in which the logical records were generated and re- 
corded, occasional adjacent pairs of records had identical 
recorded times. The hours and minutes portion of the time 
data was converted to units of seconds during the data reduc- 
mireny transeription process. 

Each logical record included an event code which corre- 
sponded to one of four specific phases within the tracking/ 
firing sequence. Event code 2032 specified that the gunner 
was in the normal tracking mode. Code 2033 reflected That 
The test director had given the gunner a command to initiate 
a firing sequence. The gunner's engagement of the fire con- 
trol lead lock switch was signified by event code 2035. The 
actual firing event corresponding to the gunner's trigger 
pull was designated by event code 2038. Each of these event 
codes was transcribed to the analysis tape without altera- 


HOD. 
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The position-location of the simulated target was re- 
ЕКС ІП Faw dara form from instrumentation connected to 
mime mirror drive pedestal. The back and forth movement 
pattern of the mirror drive resulted in raw data which 
ranged from -48 To -4048. The units of the raw data had no 
phywercal significance as initially recorded; however, after 
being transformed, the data was expressed as an angular 
measure with a range of O to 2m radians. During the conduct 
ШІ 1Пе Testing, it was noted that there was a slight amount 
of mechanical slack within the mirror drive pedestal which 
resulted in a minor deviation between the actual location of 
the mirror drive and the location which was being recorded by 
The test instrumentation. Personnel at APG in charge of the 
Test instrumentation were able to quantify and record this 
deviation thereby allowing the necessary correction of the 
mirror drive data to be accomplished during the data trans- 
mamon process. A plot of the transformed mirror drive 
position data is Shown in Appendix A, Figure 2. 

The last elements of the test data contained in the logi- 
ШЕШІ Гесога which required transformation were The gun Tube 
raw azimuth and elevation data. Gated TV video cameras were 
Mo record The dun Tube position data in units of "counts". 
ШЕНЕП 159 tne beginning of each Tracking sequence, a calibra- 
ЕИО КЕБ gated ТУ instrumentation was performed to obtain 
the proper alignment between the TV instrumentation and the 
ШОШО ЛШ ОТ The tank gun Tube. The data resulting from The 
calibration was recorded on the magnetic tape record. This 


calibration data was utilized to reexpress the azimuth and 


20 





БЕГЕ ОИ count! data as centimeter displacement relative 
To The target. IT was noted earlier that the M60AI tank 
КЕУИП У canted 90.2 mils in the counterclockwise 
econ relative To the Target. In order to establish a 
common coordinate system between the gun tube TV camera in- 
strumentation and the target, a transformation was used to 
rotate The camera coordinate system 90.2 mils clockwise. 

The final step was to convert the azimuth and elevation data 
СОСЕ 5650 in centimeters to an angular expression which 
would facilitate the planned analysis of the data. Since 
the distance between the target and the gun components was 
known and remained constant throughout the testing, The 
desired Transformation To mils was accomplished by dividing 
The centimeter deviation data by the Known distance To ob- 
tain radians. The radian measure was subsequently converted 


ШӘСПІ |5. 


С. DATA TRANSCRIPTION AND SMOOTHING 

As a consequence of the high sampling rate at which the 
test data was collected and recorded, approximately 25,500 
logical records of raw data were created for each individual 
tracking run. Because of the enormous amount of data which 
was generated during the series of tracking runs, it was nof 
Sons idered practical to analyze each Individual logical re- 
Sao data. Analysis of The data had to proceed under The 
Semomnaint imposed by computer core storage capacity, yeT 
sol be reduced in volume fo the point That significant 
name ion was lost. The amount That the data could be con- 


densed without significant information loss was largely 
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subjective and therefore did not lend itself to a unique 
БӨГЕТІ оп. 

Two separate data tapes were created using the computer 
MN previously described. Each of the tapes contained 
Со регла ипо To tracking runs 36 through 142 but they 
differed in the amount of data which each contained. The 
initial tape (tape |) created from the raw data tape con- 
tained gun azimuth and elevation tracking data and mirror 
drive position data, along with corresponding time data. 
LFogical records pertaining to times when the lead enable 
Switch was engaged and those corresponding to simulated fir- 
Ings were not transcribed. Thus, only logical records which 
feerained 10 normal tracking were transformed and copied to 
The new tape. The volume of the raw data was reduced by 
meecessing and Transcribing every 22nd logical record of the 
tracking data. Since the raw data was recorded at 106.55 
logical records per second, transcribing every 22nd record 
reduced the data to a sequence occurring approximately 4.84 
times per second. The rationale for choosing each 22nd re- 
Nas тс Obtain Tracking data points at approximately 0.2 
second intervals. A 0.2 second interval between logical re- 
cords was chosen because previous fire control testing at 
APG, similar to that being investigated in this thesis, had 
Enmplloved a data base collected at 0.2 second intervals. Sub- 
Sequent analysis of the data collected af this rate had been 
practícal in terms of the amount of computer core storage re- 
NE and also had provided a satisfactory record of the 


ces of the Tracking process. 


за 





A preliminary analysis of the test data included plotting 
gun tube azimuth and elevation data versus time for each 
us run. the plots obtained for tracking run number 6 
Emo representative of those obtained for other tracking runs. 
These plots are shown in Appendix A (Figures 3 and 4). The 
data transcribed on the initial data tape was also printed so 
une character of the data could be examined for each 
logical record. As will be explained in detail in the section 
concerning Data Analysis Methodology to follow, the relation- 
ship between varying filter delta time intervals and mean gun 
Tube azimuth error rates was analyzed. Based on the analysis 
of the data contained on tape |, it was concluded that a sub- 
sequent tape (tape 2) should be produced from the raw data 
Beer the decision to create a second tape for analysis pur- 
poses was motivated by the need to obtain greater resolution 
for the optimal time span for the tracking rate filter than 
imme O.2 second ¡interval would permit. Additionally, empirical! 
evaluation of the azimuth and elevation data revealed erratic 
fluctuations in adjacent logical records which suggested the 
need to perform a smoothing process on the raw data To 
eliminate the observed aberrations. 

Appropriate changes were made in the data reduction pro- 
gram so that the gun tube azimuth and elevation data on every 
third logical record was processed by a data smoothing routine. 
Dressing each third logical record of the raw data tape re- 
ed in a data recording density which averaged 35.52 logi- 
records тог each second of tracking Time. ln order to 


Miner. reduce the number of logical records recorded per 
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EB omg, Only each second logical record of the smoothed 
data was transformed and subsequently transcribed to magne- 
mecmudpe 2 TOr analysis. The transcription of each second 
record of the smooth data sequence resulted in a recording 
rate which averaged |7.75 logical records per second. As 
previously noted, the time between adjacent logical records 
Mas not constant and for this reason the time interval be- 
Tween subsequent logical records on the transcribed tape 
was either 0.05 seconds or 0.06 seconds. Tape 2 included 
ШЕБЕР Гус position data, and event code data in addition 
to time, and gun tube azimuth/elevation data. The entire 
Tracking, command to fire, lead insert, and firing sequence 
was included in the data transcribed to the second analysis 
tape, whereas tape | included only logical records pertain- 
ing to perids when the gunner was in the normal tracking 
разе. 

The smoothing routine utilized in the data reduction pro- 
gram to smooth the gun tube azimuth and elevation tracking 
data used the technique of taking running medians of Three 
until convergence. The method was developed by John W. Tukey. 
Ш-- ггггап coding of the smoothing algorithm was taken from 
Еа еа "Interactive Dala Analysis” by Donald R. 


Исмет [4]. 


s ANA ENS ۹ ۶۷۲ 


The purpose of This section is to describe the analysis 
which was conducted regarding the optimal "delta time" to 
ШІ а тапк Fire control tracking rate filter. Several 
different analysis approaches were utilized in this thesis. 
Each approach will be described and will be followed by the 


conciusions which were drawn from the analysis. 


٠۰٠2 ۱۱۱۷۸۷۹۷۰۲۳۴۶۸6۴ ۱۱۷ ERROR 

Tape | was used as the data base for the analysis of 
Samer mean azimuth tracking errors as a function of delta 
time. The procedure employed was to investigate how the 
Sumer azimuth tracking error changed as The delta Time 
(averaging time) was increased from zero to a maximum of 
ПЕРСО Тате у 10.25 seconds. Тһе reader will recall That 
because of the way the time data was initially recorded 
there was not a constant time interval between each logical 
record of data. The time between the logical records used 
in this analysis was either 0.20 seconds or 0.21 seconds. 
The delta time steps used in the averaging process were at 
either of these time intervals. For the sake of clarity The 
uc Kmg discussion will consider that the time between 
ШЕ ШШШ уе Logical records was a constant 0.20 seconds. 

Two different procedures were used fto calculate mean 


amet azimuth tracking error. The first procedure (pro- 
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cedure |) employed was to start at the Бе ОСИП one 
tape containing the test data and calculate the mean gunner 
БЕНШІШІШ tracking error by grouping the logical records into 
М-тиріев where N was increased from I to 50 in increments 
of 1. Considering data records in groups of N-tuples corre- 
sponded to a delta Time averaging interval of (N-1)(0.20 
Seconds), depending on the assigned value of N. As an ex- 
Eu итеп N equalled |, the tracking error from al! logical 
records was taken individually and a mean was taken of the 
entire group. This corresponded to a zero delta time. When 
N equalled 2, the mean error was computed for each adjacent 
Fogi ca | ord pair and the mean of all the pairs was taken 
٣۰۰۰٠٠٠-٠٠٣-٢٦٦٢ error for a delta Time corresponding to 0.2 
Ends. The effect of This averaging procedure was To 
eap frog" through the data by grouping logical records | 
EN Ss: a pair, 3 and 4 as the second pair until finally K-I 
ЕШ ҚК сопотітитей the final pair. An analogous "leap frog" 
effect occurred when the logical records were grouped as 
Peres on Through 50-tuples. Logic in the Fortran analy- 
Sis program detected time breaks in the normal tracking data 
caused by simulated firings and skipped these areas in the 
averaging computations. A general description of the Fortran 
Program used to accomplish the foregoing is located at 
peeendix B, flowcharts | and 2. 

The second method (procedure 2) employed to compute mean 
Fan tracking error was to calculate the error for in- 


Greasing delta time intervals starting with the last logical 





record prior to a simulated firing sequence, proceeding back- 
ШЕ ІП time. the rationale for using this scheme was based 
ФР баева thar the gunner's tracking performance was most 
ИШ DEST Chad minimum error) prior FO a simulating firing. 
ШИН cori gating mean Tracking errors in the critical time 
just prior to firing it was hoped that some relationship be- 
Den mean Gunner error and delta time could be established 
which was not detectable by the first procedure employed. 
Delta time was increased in increments of 0.2 seconds by in- 
cluding | additional logical record in the averaging process 
each time the computation was repeated. Mean errors were 
calculated for each tracking segment prior to a firing se- 
quence and then These means were combined to form an overall 
mean for the entire block of data. Flowchart 4 gives a 
description of this procedure. 

At this point it is appropriate to recall that the analy- 
sis being conducted involved searching for an optimal delta 
time to assign to the tracking rate filter in a tank fire 
control lead prediction algorithm. Lead prediction as used 
in the context of this thesis involves the problem of deter- 
mining the angular amount by which to augment the gun firing 
329-٠٠٣٦٠٦ order toWeompensate for the lateral motion of a 
femme target. The point to be made is That Тһе solution of 
the lead prediction algorithm provides a correction То azi- 
muth, not elevation. Because of this consideration the 
analysis conducted in this thesis is keyed on searching for 
an optimal delta time which will minimize azimuth firing 
ЕЙ | Errors Occurring in elevation are not considered in 


the analysis. 





A scatterplot of delta time versus azimuth error in mils 
meaeenalysis procedure | and 2 is shown in Figures 5 and 6 
K r CT е ihe mean tracking error in each of these scat- 
ТРИО decreases as delta time increases. During the track- 
ELO -gsencee The gunner’s failure to perfectly track the tar- 
ДЕТИ а= Пошта results in positive azimuth tracking errors 
Wien he OVerleads the target and negative errors when he 
underleads. With increasing delta time these positive and 
negative errors tend to cancel themselves resulting in a mean 
Tracking error which approaches zero, Comparison of the mean 
azimuth tracking errors of the two procedures for delta time 
values between zero and approximately | second shows that 
The mean errors are less for procedure 2, In procedure 2 
me errors were computed starting just prior to the firing 
sequence moving backward in time. This suggests that the 
gunner was tracking the target with less = е н 1 7٣۵۹١ 
azimuth just before the firing sequence than at other times 
in the tracking sequence. Each gunner who participated in 
ШИЕ (55! was Told prior to the initiation of testing that he 
ШОШ Га track the target as conscientiously during non-firing 
۶۰-٣٠٠٠٠٦٠٦٦٦٦ ٦١٢٥٥٠۶ as during The firing sequences. Even 
though each gunner was given these instructions, the lower 
TT SiON just prior To a firing tend to reveal fhe human 
Eon to perform most effectively during Those Times 
he perceives to be most important. 

iemeonclusion to be drawn from the foregoing analysis 
simply that increasing delta time periods increase the 


chances that the gunner's positive and negative azimuth 
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ИИТИИ errors will cancel themselves, thus resulting in a 
uude dgSeppeberending to zeros This result seems to 
argue that perhaps only large delta times should be con- 
SESA TOT USE MT the fire control lead prediction algorithm. 
This consideration will be addressed again in the analyses 


по 91 Гом. 


Б. AZIMUTH TRACKING RATE ERROR 

Ais Of The mean azimuth Tracking rate error as a 
function of delta time was accomplished in a manner very 
similar to that described for the mean azimuth tracking error. 
ЕЗИ tape | was used as the data base for the analysis. 
Analysis procedure | and 2 were modified so that mean azimuth 
Tracking rate errors were calculated where mean azimuth 
error had been previously computed. The modified versions 
of procedures | and 2 will be referred to as procedures 3 
and 4 respectively. 

As was the case for its counterpart, procedure 5 grouped 
uu auscalverecerds of data into N-tupes where N was Iterated 
uuo 90 in | unit Ітсгетепіс. Azimuth Tracking rate 
error was computed between adjacent logical records by sub- 
ШЕН Не еггог for the ¡"th record from the (i'th + 1) 
Nord error and by dividing this result by the Time between 
Mio records. Delta™fime intervals were increased accord- 
Mao the expression (N)(0.20 seconds). As an example when 
N equalled | the error rate calculation was: 


2- | n DE و‎ 07 L 
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where К ional muinmber of logical records in data base 


A 


time between adjacent logical records. 
For a delta time interval of 0.40 seconds, which would be 


the case when N equalled 2, the computation was: 


2- | 3-2 4-3 5-4 (K-1)-CK-2) (K)-(K-1) 
— — a en т 
2 
(К+ 1) Š 
2 


The general scheme shown above was repeated for delta time 
calculations through N equals 50. Whenever the program logic 
detected a time break in the tracking sequence this interval 
was skipped and the rate calculations were continued for the 
logical records following the break. A description of the 
procedure 5 error rate computation is given in flowchart 5. 

As in analysis procedure 2, procedure 4 located the last 
(раіса!) record prior to a simulated firing sequence and began 
peor calculations starting at that point proceeding backward 
in time. As before, the size of the N-tuple in a given intera- 
tion of the analysis procedure governed the length of the 
delta time period. In general the delta time period was (N) 
(0.20 seconds) for any given N. The value of N was incre- 
mented in steps of | from an initial value of | to a maximum 
EU: 95009 The rationale for beginning error rate computa- 
tions just prior to a simulated firing sequence was the same 
BR Gegen the discussion of analysis procedure 2.  Pro- 
БОО Ге 4 іс described in flowchart 5. 

The results of the mean azimuth error rate computations 
for procedure 3 and 4 are shown in Figures 7 and 8 respec- 


tively. The plots clearly indicate that the error rate 
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decreases as the interval of delta time increases. This re- 
sult is similar to that which was obtained in the case of 
The mean azimuth error. The explanation for the decreasing 
aia ga ln les in the fact that as delta time becomes 
greater, the positive and negative tracking rates begin to 
comcel Themselves resulting in a mean tracking rate error 
which approaches zero. It is interesting to compare the 
Slopes Of ihe mean error rate scatterplots (Figures 7 and 8) 
with the slopes of the mean error scatterplots (Figures 5 and 
6). The slopes of the error rate plots are greater than those 
for the mean error plots. The difference in the steepness 
of these slopes can be explained as follows. Positive and 
negative azimuth tracking errors are the result of tracking 
first in front of the target and then behind the target. In 
order to generate positive and negative azimuth tracking 
error rates it is not necessary to track back and forth 
across the target. Positive (negative) error rates can be 
caused by nonconstant movement of the turret even though 
throughout the period that the errors are generated the gun 
Tube is always leading (lagging) the target. Based on the 
foregoing it seems reasonable to expect a higher rate of 
error cancellation of mean azimuth error rates, as compared 
to those of mean azimuth errors, for a given interval of time. 
ИШ И ت۱06‎ ٣٣٢٣٣680 azimuth tracking rate errors 
are smaller for procedure 4 than for procedure 5. This sup- 
ports the earlier observation That the gunner tended То 
uu he target with less error just prior to a firing se- 


quence than during those times when a firing was not eminent. 
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КИСАЕ ШБ Of analysis procedures 35 and 4 did not reveal! 
any significant information regarding the optimal value for 
a delta time that had not been observed in analysis proce- 
dures | and 2. All four of the analyses demonstrated that 
Wan error could be reduced by increasing the length of the 


delta time interval. 


BED TALZTTARGET AZIMUTH ERROR (MODEL 1) 

The data base used for the analysis discussed in this 
section employed the data recorded on tape 2. This was the 
second tape which was transcribed from the magnetic tape 
containing the raw test data. The time interval between the 
(Са records on tape 2 was either 0.05 seconds or 0,06 
seconds. The data contained on tape 2 represented each of 
Mee ourn phases of each tracking run. This is contrasted 
Em pe i which contained only data pertaining To periods 
when the gunner was in the normal tracking sequence. 

The method of analysis described in this section was the 
5۴۰۰۰.٠٦٦ +ح)‎ +٦ quantity azimuth error as function of 
miter rate, gunner error, target rate and delta time, This 
was done in the following manner. 

Der step in the analysis process was to establish 
common coordinate system in which to express The Turret 
БИН 2206! positions. This step was necessary In order To 
be able to define relationships between target motion and 
turret motion which would permit computation of total target 
ЕПТІ error as functions of these motions. 

Among the elements of data contained in each logical re- 


cord of tape 2 was the mirror drive position expressed in 
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Ense cine range of тһе mirror drive position data was O 
to 2m radians. The mirror drive position data represented 
the position of the target (the laser spot) as it moved back 
ОШОО across the target screens Since the target motion 
was Sinusoidal, taking the sine of the mirror drive data 
resulted in transforming the position of the target to a 
scale ranging from -| to +1. The radius of the projection 
path across the target screen was computed to be 60.323 
centimeters therefore by multiplying the converted mirror 
ave position data by this value the target motion е оао > 
expressed in centimeters with a range of -60.323 to +60, 3235, 
In order to have a coordinate system with zero as an origin, 
The scale was shifted in a positive direction by adding 
S25 To each coordinate value. The result was a coordinate 
system which reflected the target motion on a scale ranging 
from 0 to 120.646 centimeters. The corresponding range of 
this coordinate system expressed in mils was 0 to 38.165 mils. 
The azimuth tracking data recorded as an element on each 
logical record indicated the deviation of the gun sight from 
пе nrarget in mils. The sign convention employed when This 
ааа иа Initially collected was fo assign a negative value 
to the deviation when the sight was to the right of the Tar- 
get and to assign a positive value when the sight was To The 
left of the target. With the knowledge of the target loca- 
tion and with the known deviation between the gun sight and 
mem iarge: it was possible to locate the position of The 
etiam the target coordinate system, Expressed in symbols 


the relationship was as follows: 
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where 
EOS dinate Of The gun sight 
Im СОО О nere of The target 


2 = deviation in azimuth between the gun sight 
and the target 


moving achieved The goal of creating a coordinate system 
common to both the sight and the target the next step in the 
analysis was to derive an expression for target azimuth error 
which would be a function of the gun and target parameters, 
as well as delta time. 

Де аке Of Clarity it is appropriate To point out 
that the gun sight rate computed in this analysis is equit- 
able to the turret rate. The only time when these quantities 
would not be identical would be for a moment of Time directly 
following the gunner's activation of the lead lock enable 
switch. That period of time is not pertinent to this analy- 
sis therefore gun sight rate and turret rate should be con- 
sidered to be synonymous. 

In this analysis the model used to compute the total 
target azimuth error was: 

F A ) UTUQGGAT)TOF 9 B = (CFITOF 


O Oj FE OJO © 


where EC = Total target azimuth error (mils) 
G(A+) = Mean turret rate delta time before 
lead insertion (mils/sec) 
TOF = Projectile time of flight (sec) 
ӨШ = nner error n az rmt Әс firing mild s 
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F 


Бал orge hate during =the TOF (mi ls/sec) 


At = Denotes that the parameter is a function 
of delta time (sec) 
The azimuth error computed using the model above will be 


positive when The target miss occurs to the right of the 
target and will be negative when the miss occurs to the left 
КОШ Оа Сет. The first and last terms in the equation will 
cancel themselves when the mean turret rate before lead in- 
Elton equals the mean target rate during the projectile 
time of flight. When these two rates are identical the lead 
computed by the fire control prediction algorithm will equal 
Шеше lead required to hit the target. The gunner azimuth 
error at the time of lead insertion was used in the model as 
a substitute for the gunner azimuth error at the time of 
Trigger pull. The actual gunner azimuth error at the time 

of firing could not be computed from the data which was trans- 
scribed on analysis tape 2. After lead insertion the gun 
wbc azimuth no longer coincided with the sight azimuth thus 
NU-SposiTyon of the sight relative to the target could not 

be determined. All of the computations performed in this 
Thesis used the gunner azimuth error at lead insertion as a 
surrogate for the actual gunner azimuth error at the time of 
pir ing. In this analysis although it was not possible to 
determine how closely the gunner azimuth error at lead inser- 
tion approximated the gunner azimuth error af firing, it seems 
КОИ ЕЗ that any error caused by This approximation is пої 
КОС tank gunners are taught in their gunnery instruction 
to obtain an estimate of the target movement rate prior to 


ШЕ Іпсегтіоп by tracking the target center of mass. After 
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lead insertion has been accomplished the gunner resumes track- 
ШІ СОШ the target center of mass until Иа Nes mo co red. 
Since the gunner is theoretically tracking the target center 
Brass ат both lead insertion and at firing, The error at 
These two distinct times should be highly positively corre- 
lated. 

The projectile time of flight used in this analysis was 
1.0699 seconds. This parameter value was obtained from 
AMSAA and it represents the expected length of time required 
for a projectile to travel a distance of [500 m. 

Flowcharts of the Fortran program used to accomplish the 
computations described in this section are located in 
ТЕЛДЕ В. flowcharts 8 through I5. The actual azimuth 
ENocpecomputation is accomplished in the subroutine described 
КУ lowchart 12. 

iNET DOINGS Of The total target azimuth error versus 
delta time is shown in Figures 9 and 10. Figure 9 shows the 
mean azimuth error over a delta time interval ranging from 
0.06 seconds to 5.59 seconds. The delta time interval was 
expanded from its lower value То its upper value by including 
Иоса епа logical record in the error calculation each 
time the computational loop was repeated. This resulted in 
an incremental delta time step of either 0.05 seconds or 0.06 
seconds. The maximum delta time value of 5.59 seconds repre- 
sented the inclusion of 100 logical records in the azimuth 
EN Uai on. The plot shown in Figure IO reflects the 
azimuth error result when a maximum number of 55 logical re- 


cords were used in the computation. Limiting the maximum 
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number of logical records (used in the computation) to 33 
ШОШ е effect of magnifying The detail of the scatterplot 
xs data points occurring in the delta time interval 
ШЕШІП 9,05 Seconds to 1.81 seconds. 

ШЕНШЕе 9 51оМ5 that The computed azimuth error is a mini- 
MS e delta time equal to approximately 0.68 seconds, but 
Eer esses for all delta times which are greater than 0.68 
Emus The point where the minimum error occurs is shown 
Mires mer detail in Figure |0. From this figure the delta 
time corresponding to minimum target azimuth error can be 
determined to be equal to approximately 0.70 seconds. 

In order to gain a clearer insight concerning why the 
azimuth error minimum occurred at 0.70 seconds a comparison 
was made of the mean deviation between the turret rate before 
lead insertion (as a function of delta time) and the target 
rate during the projectile time of flight. It can be seen 
from the target azimuth error equation that azimuth error is 
sensitive to differences between these two rates. Figures 
regel ZZzsnon The results of the rate comparison. Both of 
These figures show That the minimum deviation between The 
TT Tale and target rate occurred at a delta time inter- 
Sar 0770 seconds. Comparison of Figure |! with Figure 
9 shows that the two plots are essentially identical with the 
КЕК ЫЕ enethat minimum ordinate of Figure ll is slightly less 
than that of Figure 9. This difference is caused by The in- 
Sion ot the gunner azimuth error at lead insertion In The 
target error computation whereas this error was not included 


Mine rate comparison. 
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The results of this analysis demonstrated that azimuth 
error is highly sensitive to the error which results because 
uere: rate before lead insertion: is not equal to the 
I Target rate after firing. iunc ona steh eS 
Sieamtiicd above, the error resulting from the gunner's 
Tailure to have the sight reticle on the target in azimuth 
ЕШ Ее rime от firing was not great. This suggests that the 
magnitude of the target azimuth error is very dependent upon 
umuecerrecTiveness of the fire control lead prediction algo- 
ТІ him. 

After the analysis described above had been accomplished, 
a critical appraisal of the target azimuth error model was 
made for the purpose of identifying improvements which could 
be made in the model structure. Keeping in mind that the 
primary purpose of the entire analysis was to investigate 
ШЕ principal dynamics Of azimuth error as a function of 
delta time, it was considered desirable То redefine The model 
in Terms more descriptive of the real world. The redefini- 
tion of the model and the subsequent analysis on The eluci- 


Hed" model are described in The next section. 


ЕР TOTAL TARGET AZIMUTH ERROR (MODEL 2) 

mW analysis was Similar in many respects to the anafy- 
Sis described in the preceding section. The analysis was 
ОООО ОТ а on the data recorded in tape 2 and made use of The 
E ccgures required To construct a target coordinate system 


based on the mirror drive position data. 
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The model used to compute total target azimuth error as 
a function of deita time, target rate, turret rate and gunner 
error was a summation of three sources of azimuth error. The 
first source of error was that error caused by an incorrect 
lead computation resulting because the turret rate was not 
БОША o rne target гате in a variable delta time interval 
prior to lead insertion. The difference between these two 
ЕОС ир Тес Бу the time of flight of the projectile is 
o he error due to incorrect lead. This first error 
Lx cs dhead prediction error. The second source of error 
identified in the model is the error created as a result of 
E ssuuSE' s failure to place the sight reticle on the target 
meet Ot Mass in azimuth at the time of firing. This error 
will be referred to as gunner azimuth error. The third and 
final source of azimuth error included in the model is that 
error caused by a target rate after firing which is different 
from that computed during a given delta time interval prior 
to the gunner's insertion of lead. The difference between 
These Two target azimuth rates multiplied by the projectile 
time of flight will be called the target induced error. The 
target induced error reflects the error caused by evasive 
movement of the target after firing which was not predicted 
based on the target's movement in a delta time interval prior 
to lead insertion. AS an example, it would be logical Tto 
ns that target induced error would be high for target 
Mehicles possessing high mobility/agility. The model com- 


EET GT the errors just described can be expressed as follows: 
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EREMO = АГАТ) + B * C(AT) 
DR = GCAT))TOF + B + CE =- DCAT)) TOF 
EE UT TOF Fr B 
where 
HS oral target azimuth error (mils) 
oni) = bead prediction error (mils) 
rss jn (error In azimuth at firing (mils) 
С(А+) = Target induced error (mils) 


D(At) = Mean target rate delta time before lead insertion 
(mils/sec) 


ЕШ Mean turret rate delta time before lead insertion 
anms sec) 


lU = rrojecrile time of flight (sec) 
ПОЕ ИСЗЛЕ 5 гает rate during the TOF (mils/sec) 


AT = Denotes that the parameter is a function of delta 
time (sec) 


The target azimuth error expressed by this model will be 
mesitive if the target miss occurs to The left of the target 
and will be positive for misses to the right of the target. 

From an examination of The model above it can be seen that 
mies product D(At)TOF occurs twice and thet it cancels itself. 
After the cancellation of these terms the resulting model 
is identical to the model used in the previous analysis of 
ШЕ Го та aget azimuth error. The value of expressing The 
ШЕ | ІТ is done in this analysis Is to permit Тһе total 
Peden azimuth error to be segmented into three distinct 
КООШО СТ error. By quantifying The azimuth error in Terms 
Oris unique sources it was hoped that insights could be 
Samed concerning which of these three sources contributed 


Seemesignitficantly to the total error. 
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General flowcharts of the Fortran program used to imple- 
ment the azimuth error model are located in Appendix B, flow- 
TIS $ Through |5. 

The results of the analysis described in this section 
be presented by using scatterplots that display Total 
meget azimuth error sources as a function of Increasing del- 
ta time, Because the analysis objective was to gain insights 
regarding the relationship between the target azimuth error 
and delta time, emphasis was placed on the general trend of 
The target error as delta time increased rather than on the 
magnitude of the error. 

The scatterplot of target induced error versus increasing 
delta time is shown in Figures 13 and 14. These two figures 
were created from the same data; however, the plot shown in 
Figure 14 reflects a delta time period which is approximate- 
K í jd of That shown in Figure 13. Figure 13 shows a 
general positive trend in the target induced error beginning 
at a delta time of approximately 1.0 seconds. This positive 
trend results because as delta time is increased the mean 
traverse rate of the turret becomes increasingly less repre- 
Sentative of the target rate after firing. AS was observed 
I Tracking rate analyses (Section VB), the mean azimuth 
Tracking rate tends to converge to zero with increasing delta 
NEC Xgure 14 shows that the target induced error global 
minimum occurs at a delta time equal to approximately 0.95 
Eus This indicates that for the data used in this 
meh sis the target induced error can be minimized by select- 


ШЕШ Ше | Та time of 0.93 seconds. Because fhe azimuth error 
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modet includes two other error sources in addition to target 
ІШШСеЭей еггог, іт cannot be concluded that a delta time of 
uu sdurarromn minimizes total azimuth target error. 

A scatterplot of the mean deviation between the turret 
Ba iq the Target rate prior to lead insertion is shown in 
Figure 15. The deviation in the rate is plotted versus in- 
астас авта time. The deviation multiplied by the pro- 
КЕСИБ fime Of Tlight (a constant) represents the source of 
error earlier identified as lead prediction error. The de- 
Smedsing trend of this plot reflects an inverse relationship 
between increasing delta time and the rate deviation. As 
the delta time interval is increased the instantaneous track- 
ing errors made by the gunner tend to cancel themselves with 
The result that the average turret rate becomes an increas- 
ingly-better approximation to the actual target rate. 

The two sources of target azimuth error discussed above 
represent opposite error trends as delta time is increased. 
The lead prediction error decreases with larger delta time 
whereas the target induced error increases. The results of 
the summation of these two sources of error are shown in 
ОЕ eS бапа 17. This scatterplot shows that the Target 
induced error clearly dominates the lead prediction error. 
BEEN c0 nePigure 17 the global minimum of the sum of These 
errors occurs at the delta time equal to approximately 0.70 
seconds. 

mI Pd source of target error not considered in Fig- 


Messe and 17 is The gunner azimuth error at firing. This 
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error is independent of delta time. For the reasons explained 
previously, gunner error at lead insertion was used as a sub- 
stitute for the actual gunner error at firing. When the 
three components of the target azimuth error model are summed 
and plotted versus delta time, the results are as shown in 
mares |G and 19. By comparing these figures with those 
shown іп 16 апа 17 it is apparent that gunner azimuth error 
Aina C as modeled by the gunner error at lead insertion) 
ИШЕ major contributor to the total target azimuth error. 
The plots shown in Figures I8 and I9 are identical to those 
shown in Figures 9 and 10 thus verifying that the models 
used in analysis method | and method 2 were identical except 
for the manner in which they were expressed. Both models 
Ear a delta Time corresponding to approximately 0.70 
seconds results in the minimum total target azimuth error. 
Me ree sources of error included. in This Target 
error model, the error term which contributes most to the 
ial Ta GCE azimuth error is the target induced error, The 
scatterplot shown in Figure 20 represents the mean azimuth 
error minus the target induced error. Comparison of this 
ШЕШ е УГЕ Figure |8 shows That the Target induced error 
clearly dominates the other factors in the model. Without 
the target induced error in the model the total target azimuth 
error is drastically reduced. Considering the dominance of 
the target induced error over the other error sources, it is 
apparent that total target azimuth error can most effective- 


ly be minimized by reducing the target induced error. 
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VI. CONCLUSIONS 


In this Thesis several analysis techniques were developed 
for the purpose of analyzing the tank fire control test data. 
The analysis objective was to investigate the optimal time 
(delta time) to assign to the tracking rate filter employed 
in a tank fire control lead prediction algorithm. An addi- 
Mona! objective was to gain insights concerning the dynamics 
of target error. This involved the identification of sources 
of target error and the investigation of the relationship be- 
Tween these sources and delta time. 

The relationship between azimuth tracking error and delta 
Time was analyzed in Section VA. Both analysis procedures | 
and 2 demonstrated that the gunner mean azimuth tracking 
error decreased as delta time was increased. The gunner mean 
azimuth tracking error was shown to be smaller just prior to 
firing sequence than at other times in the tracking sequence. 
This suggested that the gunner was motivated to track The 
Pegetemore conscientiously during those fimes just before 
a firing sequence than at those times when firing was not 
imminent. 

Analysis procedures 5 and 4 in Section VB investigated 
sumer mean azimuth Tracking rate error as a function of 
Nue. The results showed that fhe mean azimuth frack- 
ing rate error decreased as a function of increasing delta 


time, 
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The results obtained from analysis procedures | through 
4 demonstrated that the mean azimuth tracking error and mean 
azimuth tracking rate error could be minimized by assigning 
Esugecvalblues to the tracking rate filter. The insights 
gained serve to motivate the subsequent analyses which were 
conducted but did not provide any findings that isolated a 
unique value to assign to the tracking rate filter. 

Based on the initial analysis it became apparent that 
СЕЕ есине analyses should "include the portions of the test 
data pertaining to the firing sequences as well as the data 
representing the non-firing sequences. In addition it was 
concluded that it would not be possible to arrive at any 
meaningful findings relative to an optimal delta time value 
unless the interval of time between adjacent data records 
was reduced by including more records per unit time in the 
analysis. Based on these findings a second analysis tape 
was transcribed from the raw test data which included The 
additional data that was not included in the initial analysis. 

Model | described in Section VC was formulated to quanti- 
Em code varget azimuth error as a function of delta time. 
The results of this analysis procedure showed that minimum 
meget cemror occurred at a delta time equal To 0.70 seconds. 
It was determined that error produced by the gunner's failure 
mee have the sight reticle on the center of mass in azimuth 
was small in comparison to the error caused by movement of 
the target after firing that was not predicted by the lead 


wee ion algorithm prior to firing. 
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pie sic) Model 2 quantitied total target azimuth error 
Boe sum Of three individual sources of error. The three 
٠۰۰٠٠٠٠٦٦۷6 ۴٣۴ lead prediction error, gunner azimuth error, 
mudEarget induced error. The lead prediction error and the 
target induced error were investigated individually to deter- 
Be K lc ir Characteristics as a function of delta time. 

The results of the analysis employing Model 2 revealed 
an inverse relationship between the lead prediction error and 
increasing delta time. Target induced error was shown to in- 
crease as delta Time became greater. The gunner azimuth 
error was not a major component of azimuth error when com- 
pared to the lead prediction error and the target induced 
٣٣ٍ rhe target induced error component dominated The 
other two components in terms of magnitude and therefore was 
usur cendriburor To the total target azimuth error. 

Both Models | and 2 demonstrated that total target azimuth 
error was minimized at a delta time equal to approximately 
0.70 seconds. 

The value of the analysis procedures described in this 
mies is does not lie entirely in the determination That The 
miminumetotal target azimuth error occurred at a delta Time 
EIN Omceconds. This finding is no doubt dependent on The 
character of the test scenario and therefore would be expec- 
ted to change somewhat if other scenarios were used. The 
rq c ot the analysis is that Iit lends insights concern- 
ing which components of the total target azimuth error are 
most dominate and how components vary as a function of delta 


Time. 
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F hp + LGURES 


The figures depicted in this Appendix show plots of the 
raw test data and results of the analyses conducted in the 
thesis. 

Figures 2 through 4 are plots of the raw test data. 
Figures 5 through 8 show the results of analysis procedures 
| through 4 respectively. 

The results obtained from the analysis employing Model | 
are shown in Figures 9 through 12. Figures 15 Through 20 
depict the results of the analysis in which Model 2 was 
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TNE fiowcharts contained in this Appendix provide a 
verbal description of the two Fortran computer programs 
used in the analysis of the fire contro! test data. The 
analysis concerned with the computation of azimuth track- 
Ш С еггог and azimuth tracking rate error is described 
EN cweharTrs | through 7. Analysis procedures | Through 
4 are explained in flowcharts 2, 4, 3, and 5 respectively. 

ШЕ о лг е 8 Through 15 describe The computer program 


written to execute analysis Models 1 апа 2. 
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FLOWCHART |. MAIN PROGRAM 
(Analysis Procedures | through 4) 


Read data from magnetic tape |. Each logical record read contains 
the following information: 


Azimuth tracking data (mils) 

Elevation tracking data (mils) 

Mirror drive position data (radians) 

Recording time of the logical record (0.2 sec) 


Y 
Record the number of logical records which were read (N). 
Y 
specify the number of iterations (ITIMES) which are to be accom- 
plished within each subroutine. This specification will cause mean 


azimuth error and mean azimuth rate error to be computed by averaging 
over delta times ranging from approximately 0.2 seconds to the product 
of (ITIMES) (0.2 seconds) in increments of 0.2 seconds. 
| + 
Cal! Subroutine BLKERR. This subroutine computes mean azimuth 
tracking error beginning with tracking data following a trigger pull, 
working forward in time. 
Y 
Call Subroutine AVRATE. This subroutine computes mean azimuth 
tracking rate error beginning with tracking data following a trigger 
pull, working forward in time. 
y 
Call Subroutine ETUPLE. This subroutine computes mean azimuth 
Pedemingectror beginning with Tracking data just prior to lead insert, 
working backward in time. 
Y 
Call Subroutine RTUPLE. This Subroutine computes mean azimuth 
ШЕЛ КОЛО rate error beginning with Tracking дата just prior To lead 
insert, working backward in time. 
Y 
Call Subroutine OUTPUT. This subroutine provides a printed output 
of the data calculated in previously called subroutines. 
Y 
Cal! Subroutine PLOTS. This subroutine creates scatterplots of 
the data calculated in previously called subroutines. 


fid 
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FLOWCHART 2. SUBROUTINE BLKERR 


Enter loop 20 which will be iterated ITIMES starting with L=I. In- 
dex L is incremented in steps of |. The product (L)(0.2 seconds) equals 
The delta time increment under consideration for any given iteration of 
сор 20. 

Y 

Enter loop 10 which will be iterated N times as specified from the 
main program. The constant N corresponds to the number of logical records 
To be processed. 

Y 

Identify the beginning of each tracking sequence which follows a 
simulated trigger pull! by locating successive logical records which are 
separated in time by more than (L)(0.25 seconds). 


Y 
Enter loop 6 which will be iterated L times as specified by the index 
лоор 20, 
Y 


Calculate the cumulative sum of azimuth tracking error from each 
successive logical record proceeding forward in time. Maintain a count 
of the number of sums taken. 

+ 
Ех! т К оор 6 
Y 

Calculate the mean of the cumulative sum recorded in loop 6. Calcu- 

late the absolute value cumulative sum of the mean computed in this step. 
Y 
Calculate the time interval over which the cumulative tracking error 
computed in loop 6 was taken. Maintain a cumulative sum of these times. 
+ 
Ех ) +40 0 I0 
Y 

Calculate the mean azimuth tracking error from the cumulative sum 
computed in loop 10. Calculate the mean time interval over which the 
tracking error was computed. 


Y 
Record these values in vectors whose index value is equal to L. 
Y 
۹ 0 
Y 
Return 
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FLOWCHART 5. SUBROUTINE AVRATE 


Enmer eon ZO which will be iterated [TIMES starting with index L=l. 
Index L is incremented in steps of |. The product (L)(0.2 seconds) 
equals the delta time increment under consideration for any given itera- 
moneo loop 20, 

Y 

Enter loop 10 which will be iterated N times as specified from the 
main program. The constant N corresponds to the number of logical records 
То be processed. 

+ 

Identify the beginning of each tracking sequence which follows a 
Simulated trigger pull by locating successive logical records which are 
Separated in time by more than (L) (0.25 seconds). 

Y 

Enter loop 6 which will be iterated L times as specified by the index 
Soop 20. 

Y 

Calculate the cumulative sum of azimuth tracking rate error from 
each successive logical record proceeding forward in time. Rates are 
calculated by subtracting the i'th azimuth error from the j'th azimuth 
error and dividing this sum by the time between these successive errors. 
Maintain a count of the number of sums taken. 

+ 
Exit loop 6 
Y 

Calculate the mean of the cumulative sum recorded in loop 6. Calcu- 
late the absolute value cumulative sum of The mean computed in this 
step. 

+ 

Calculate the time interval over which the cumulative tracking 
rate error computed in loop 6 was taken. Maintain a cumulative sum of 
these times. 

+ 
Exit loop | 0 
Y 

Calculate the mean azimuth tracking rate error from The cumulative 
sum computed in loop lO. Calculate the mean time interval over which 
the tracking rate error was computed. Record these values in vectors 
whose index value is equal To L. 


Y 
Exit Сор 0 
+ 
Return 
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FLOWCHART 4. SUBROUTINE ETUPLE 


Enter loop 10 which will be iterated N times as specified from the 
main program, 
+ 
Identify The last logical record in each individual tracking sequence 
by locating successive logical records which are separated in time by 
more than 0.25 seconds. Create a vector containing the location number 
of the logical records identified. 
Y 
Enter loop 30 which will be iterated ITIMES starting with index L=l. 
Index L is incremented in steps of |. The product (L)(0.2 seconds) 
equals the delta time increment under consideration for any given itera- 
tion of loop 30. 
+ 
Enter loop 20 which will be iterated IA times. IA corresponds to 
The number of individual tracking sequences previously identified. 
+ 


Encounter logic to determine when the end of an individual tracking 
Sequence is reached. When the end of a tracking sequence is identified, 
transfer back to the top of loop 20. 

Y 

Enter loop 25 which will be iterated L times as specified by the in- 
uxor loop 30. 

Y 

Calculate The cumulative sum of azimuth tracking error from each 
successive logical record proceeding backward in time from the last 
record in each individual tracking sequence. 


Calculate the cumulative sum of the time between each successive logi- 
cal record wherein tracking data was summed in the preceding step. 
Y 
2۶۱۰۱۰۹۱۲۷۹735 
+ 
Calculate the mean of the cumulative sums recorded in loop 25. Cal- 
culate the absolute value cumulative sums of the means computed in this 
step. 
y 
Ehle OOS ZU 
+ 
Calculate the mean azimuth tracking error from the cumulative sum 
computed in loop 25. Calculate the mean time interval over which the 
tracking error was computed. Record these values in vectors whose index 
value is equal to L. 
+ 
Кетигп 
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FLOWCHART 5. SUBROUTINE RTUPLE 


Enter loop 30 which will be iterated ITIMES starting with index L-l. 
Index L is incremented in steps of |. The product (L)(0.2 seconds) 
equals the delta time increment under consideration for any given itera- 
Tien or loop 20% 

Y 

Enter loop 20 which will be iterated IA times. 1А corresponds to 
The number of individual tracking sequences previously identified in 
Eubrourine ETUPLE. 

+ 

Encounter logic to determine when the end of an individual tracking 
sequence is reached. When the end of a tracking sequence is identified, 
transfer back to the top of loop 20. 

+ 

Enter loop 25 which will be iterated L times as specified by the in- 
dex of loop 30. 

Y 

Calculate the cumulative sum of azimuth tracking rate error from 
each successive logical record proceeding backward in time. Rates are 
calculated by subtracting the i'th azimuth error from the j'th azimuth 
error and dividing this sum by the time between these successive errors. 
Maintain a count of the number of sums taken. 

+ 

Calculate the cumulative sum of the time between each successive 
logical record wherein tracking rates were computed and summed in the 
preceding step. 

Y 
6 6-0 6 
y 

Calculate the mean of the cumulative sums recorded in loop 25. Cal- 
culate the absolute value cumulative sums of the means computed in This 
step. 

y 
Exit loop 20 
y 

Calculate the mean azimuth tracking rate error from the cumulative 
sum computed in loop 20. Calculate the mean time ¡interval over which 
the tracking rate errors were computed. 


Y 
Record these values in vectors whose index value is equal То L. 
+ 
Exi T 50 
+ 
Return 
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FLOWCHART 6. SUBROUTINE OUTPUT 


Create a printed output of mean azimuth tracking 

corresponding delta time data computed in Subroutine 
Y 

Create a printed output of mean azimuth tracking 

corresponding delta time data computed in Subroutine 
Y 

Create a printed output of mean azimuth tracking 

corresponding delta time data computed in Subroutine 
Y 

Create a printed output of mean azimuth tracking 


errors and the 
BLKERR. 


rate errors and the 
AVRATE. 


errors and the 
ETUPLE. 


rate errors and 


the corresponding delta time data computed in Subroutine RTUPLE. 


y 
Return 
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FLOWCHART 7. SUBROUTINE PLOT 


Create a scatterplot of delta time versus mean azimuth tracking error 
as computed in Subroutine BLKERR. 
Y 
Create a scatterplot of delta time versus mean azimuth tracking rate 
error as computed in Subroutine AVRATE. 
Y 
Create a scatterplot of delta time versus mean azimuth tracking error 
as computed in Subroutine ETUPLE. 
+ 
Create a scatterplot of delta time versus mean azimuth tracking rate 
error aS computed in Subroutine RTUPLE. 
+ 
Return 
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FLOWCHART 8. MAIN PROGRAM 
(Analysis Models | and 2) 


Read data from magnetic tape 2. Each logical record read contains 
the following information: 


Azimuth tracking data (mils) 
Elevation tracking data (mils) 
Mirror drive position data (radians) 
Event code designation 
Recording time of the logical record (0.06 sec) 
Y 
Record the number of logical records which were read (N). 
Y 
Subtract a constant K from each logical record recording time so that 
the resulting times can be represented more accurately in the binary 
number system. | 
+ 
Call Subroutine LOCPUL. This subroutine locates the first logical 
record in each firing sequence and creates a vector KTPULL (1) which 
records the location of each record identified. A count of the number 
of firing sequences detected (NPULL) is computed. 
Y 
specify the value of the constant NDELTA. This value determines The 
length of the delta time interval to be considered in the analysis. The 
delta time interval will equal (NDELTA) (0.06 seconds). 


Y 
Enter loop 5 which will be iterated once for each firing sequence of 
which there are NPULL. 
Y 


Call Subroutine SILEAD. This subroutine performs computations rela- 
tive to the location and movement of the tank turret (sight) and simu- 
lated target a variable time interval prior to the gunner's insertion of 
azimuth lead. 

Y 

Call Subroutine LASER. This subroutine performs computations rela- 
tive to the location and movement of the simulated target after the 
gunner's insertion of azimuth lead. 

Y 

Call Subroutine ERROR. This subroutine computes total target azimuth 
error as a function of delta time based on input data from subroutines 
Seen aad LASER. 

Y 











+ 
Call Subroutine ABSDIF. This subroutine computes the mean deviation 
between the turret rate and the target rate as a function of delta time. 
+ 
Enter loop 10 which will be iterated NDELTA times in increments of 1. 
+ 
Calculate the cumulative sum of the times between successive logical 
records and maintain this result in the vector STIME. 
Y 
Calculate the cumulative sum of the target azimuth errors computed in 
subroutine ERROR. 


+ 

290 ۵ 
+ 

Exide Hoop. > 
+ 


Calculate the mean time interval between successive logical records 
by using The cumulative time data in vector STIME. The times computed 
in This step represent the mean intervals in which delta time is incre- 
mented as if is increased from O to (NDELTA)(0.06 seconds). 

+ 

Call Subroutine OUTPUT. This subroutine provides a printed output 

of selected data previously computed. 


Cal! Subroutine PLOTS. This subroutine creates scatterplots of 
selected data previously computed. 
+ 
END 
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FLOWCHART 9. SUBROUTINE LOCPUL 


Enter loop 5 which will be iterated N times as specified from the 
main program. The constant N corresponds to the number of logical 
records To be processed. 

Y 

Check The event codes of each logical record to determine if They 
are code 2038. Event code 2038 corresponds to a firing sequence trigger 
pull. When a trigger pull event is detected jump out of loop 5 and 
record the logical record number where the trigger pull occurred in 
meeror KIPULL(I). 

Y 

Return to loop 5 and start the search for subsequent trigger pull 

event codes beginning at the point where the loop was exited. 
Y 

Exit loop 5 when the last logical record to be processed has been 
reached. 

Y 
RETURN 
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FLOWCHART 10. SUBROUTINE SILEAD 


Enter loop 5 for the purposes of locating the logica! record at which 
The gunner engaged the lead lock enable switch thereby inserting lead into 
The firing azimuth. Exit loop 5 with the number of the logical record 
where lead was inserted. 
+ 
Compute the elapsed time between lead insertion and the gunner trigger 
pu. 
+ 
Enter loop 15 which will be iterated NDELTA times in increments of |. 
+ 
Create a vector SPOTML which represents the coordinate of the simulated 
target in mils beginning at the time of lead insert to the time equal to 
Time of lead insert minus (NDELTA)(0.06 seconds). Hereafter, for the pur- 
pose of brevity, this time interval will be referred to as time interval 
ALPHA. 
+ 
Create a vector APRIME which represents the coordinate of the turret 
(sight) in mils in the target coordinate system for time interval ALPHA. 
Both vector SPOTML and APRIME have a maximum length of NDELTA and contain 
coordinate data for each logical record within time interval ALPHA. 
+ 
Enter loop |9 where the gunner azimuth error is computed for time in- 
terval ALPHA. Gunner azimuth error is positive if the gunner is Tracking 
in front of the target and is negative if the gunner is lagging the target. 
+ 
Ex 1 ۲۷ 
+ 
Enter loop 30. Within this loop the vector AZRTUP and SPOTRT are com- 
puted. Vector AZRTUP represents The azimuth traverse rate of the turret 
in mils/sec for time interval ALPHA. The vector SPOTRT represents the 
azimuth rate of the target in mils/sec over time interval ALPHA. 
Y 
Calculate the cumulative sum of the time interval between successive 
logical records within time interval ALPHA. Store this result in vector 
BELTMI. 


+ 
Ех ТЕ Ор O 
+ 
Enter loop 40 which is iterated NDELTA times in increments of 1. 
+ 


Calculate the cumulative sum of the deviation between the turret rate 
and the target rate in time interval ALPHA. 


Y 


QIE 





Y 
Exit loop 40 
Y 
If the firing sequence is the last one in the series of firing se- 
quences enter loop 45 and compute the mean deviation between the turret 
rate and target rate in time interval ALPHA. 
Y 
Exit loop 45 
Y 
If The variable JRITE is set equal to | enter loop 55 which will 
cause data cards with gunner target error and target apparent velocity 
to be punched. l 
+ 
|f the variable JPLOT is set equal to | create a scatterplot of tar- 
get apparent velocity versus gunner target error. 
Р 
RETURN 
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FLOWCHART 11. SUBROUTINE LASER 


Enter loop |. Compute a vector TGTLD which represents the coordinate 
of the target in mils beginning at the time of lead insertion through the 
time at which the gunner initiates a trigger pull. Vector TGTLD will be 
of variable length which is dependent on the time interval between the 
insertion of lead and the initiation of a trigger pull. Coordinate data 
is calculated for each logical record within this time interval. 

Y 
Exit loop | 
Y 

Enter loop 2 and compute the azimuth rate of the target in mils/sec 
from the time that lead is inserted through the time of trigger pull. 
store the cumulative sum of the rates computed. 

Y 
8۰ + 2 
Y 

Compute The mean azimuth rate of the target from the time of lead in- 

Seri through the time of trigger pull. 
+ 


Enter loop 5. This loop will be iterated 25 times. Compute a vector 
TGTMIL which represents the coordinate of the target in mils beginning at 
The Time of Trigger pull through the next successive 25 logical records. 
The computation of 25 target coordinates will permit subsequent computa- 
tion of the target rate through a simulated projectile time of flight of 
|.Q669 seconds. 

Y 
EXT ۶ی‎ 2 
Y 

Enter loop !0 and compute the target rate in mils/sec over an inter- 
val of time corresponding to a simulated projectile time of flight of 
1.0669 seconds. Maintain a cumulative sum of the computed target track- 
img. rates. 

4 
Ex il OCD akg 
Y 

Compute the mean target rate from gunner trigger pull through a simu- 

lated projectile time of flight of 1.0669 seconds. 
4 

Enter loop 18 and compute the target induced error from the time of 
lead insertion through the time of simulated projectile impact.  Projec- 
tile impact is assumed to have occurred 1.0669 seconds after firing. 

| 
3+8 ІСӨР ПӘ 


Y 








4 
Enter loop 20 which will be iterated NDELTA times in increments of |. 
+ 
Using 3 different computational formulas compute target error in mils 
and store the results of each computation in a separate vector. Calculate 
the cumulative sum of the target errors for each firing sequence. 
Y 
Ех 10600220 
4 
|f the firing sequence is the last one in the series of firing se- 
quences enter loop 25 and calculate the mean target error for each of the 
5 target error sums accumulated in loop 20. 


y 
۷۷۹۶ 
4 
RETURN 
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FLOWCHART 12. SUBROUTINE ERROR 


Compute the variable TARGET which is equal to the product of the 
target rate after firing times the projectile time of flight. 
y 
Enter loop 5 which will be iterated NDELTA times as specified from 
The main program. 
Y 
Compute the vector SIGHT which is equal to the product of the turret 
rate prior to lead insertion times the projectile time of flight. The 
vector SIGHT has length NDELTA. 
Y 
Compute tota! target azimuth error in mils by subtracting the quan- 
tity TARGET and the gunner azimuth error at the time of lead insertion 
from the vector SIGHT. The resulting vector is total target azimuth 
error in mils as a function of increasing delta Time. 
Y 
Exit loop 5 
Y 
If the variable IPRT is set equal to | output selected computational 
results from subroutines SILEAD, LASER and ERROR, 
y 
RETURN 
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FLOWCHART 15. SUBROUTINE ABSDIF 


Enter loop 5 which will be iterated NDELTA times in increments of |. 
Y 
Compute the cumulative absolute value difference between the turret 
rate before gunner insertion of lead and the rate of the target after 
firing. Maintain a count of the number of sums taken. 
Y 
2 + +۷۷) 
Y 
If this firing sequence is not the last in the series of firing se- 
quences return to the main program. 


Enter loop 10 which will be iterated NDELTA times in increments of |. 

Y 
Compute the mean deviation between the turret rate prior to lead 

insertion and the target rate after firing as a function of delta Time. 
Y 
Exit loop 10 
Y 
RETURN 
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FLOWCHART 14. SUBROUTINE OUTPUT 


Create a printed output of the logical record location 
lated firings. 
Y 
Create a printed output of selected variables computed 
ШЕЛІ), LASER and ERROR, 
+ 
Create a printed output of selected variables computed 
program. 
y 
Create a printed output of selected variables computed 
SIP AD. 
Y 
RETURN 
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all simu- 


Subroutines 


the MAIN 


subroutine 





FLOWCHART 15. SUBROUTINE PLOT 


Create a scatterplot of delta time versus mean total target azimuth 
error as computed tn Subroutine ERROR. 


M 


Create a scatterplot of delta Time versus the mean deviation between 
the turret rate prior to lead insertion and the target rate after firing 


as computed 


Create 
the turret 
Subroutine 


Create 


in Subroutine ABSDIF. 
y 
a scatterplot of delta time versus 
rate and target rate prior to lead 
Eh 
y 
a scatterplot of delta time versus 


error as computed in Subroutine LASER. 


Create 
nurh error 


Create 
muth error 


Create 
muth error 


y 
a scatterplot of delta time versus 
(TGTER2) as computed in Subroutine 
y 
a scatterplot of delta time versus 
(TGTER3) as computed in Subroutine 
y 
a scatterplot of delta time versus 
(TGTER5) as computed in Subroutine 
Y 
RETURN 
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the mean deviation between 
insertion as computed in 


the mean target induced 


the mean total target azi- 
LASER. 


the mean total target azi- 
LASER. 


the mean total target azi- 
LASER. 
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